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Intrinsic transverse parton momenta pr play an important role in the understanding of az- 
imuthal/spin asymmetries in semi-inclusive deep-inelastic scattering (SIDIS) and the Drell-Yan 
process (DY). We review and update what is presently known about pr from these processes. 
In particular, we address the question to which extent data support the popular Gauss model for 
the pr-distributions. We find that the Gauss model works very well, and observe that the intrin- 
sic transverse momenta in SIDIS and DY are compatible, which is a support for the factorization 
approach. As a byproduct we recover a simple but practical way of taking into account the energy 
dependence of py-distributions. 
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I. INTRODUCTION 



Intrinsic transverse parton momenta in hadrons can be probed in deeply inelastic reactions, such as SIDIS and DY, 
when adequate transverse momenta in the final state are measured. Here "transverse" means with respect to the hard 
momentum flow in the process, e.g. in SIDIS transverse momenta of produced hadrons with respect to the virtual 
photon. Transverse parton momenta are described in terms of transverse momentum dependent parton distribution 
functions (TMDs) or analogously generalized fragmentation functions [H-Q- 

Although the concept of transverse parton momenta dates back to early days of QCD [l|, 0, the field has 

received continuous interest from theory [T^-HSl and important steps in the understanding of TMDs within QCD were 
taken only recently [23T - |27| ], In particular, factorization theorems have been formulated, which extend previous work 
and ensure that the productions of hadrons in SIDIS or dileptons in DY with small transverse momenta compared 
to the hard scale factorize in hard scattering parts and universal non-perturbative objects: TMDs, fragmentations 
functions, and soft factors [26l.[27j. 

The appealing perspe ctive of the TMD approach is that it may explain single spin or azimuthal asymmetries in 



various reactions [28H30J, especially in SIDIS [3ll-l49l . DY [50l - l57j . or hadron production in e*e~-annihilations (58l - [60| . 
In phenomenological studies of such data [6ll484| one often works in lowest order QCD ("tree level") and neglects 
higher order QCD corrections, which includes soft factors as well as higher order corrections in the hard scattering. In 
some sense, the information on soft factors can be thought of as "reshuffled" into an effective description of involved 
TMDs or fragmentation functions. 

As in such studies typically many novel TMDs are involved [85T4881 ]. it is popular to assume the so-called Gauss 
model. For example, in the case of the unpolarized TMD or fragmentation function one assumes 



f?{x,p T ) = f1{x) 
D?(z,K T ) = Df(z) 



exp(-p|/{p|)) 



ir{p 2 T ) 



expi-KZ/iK*)) 



n(K 2 ) 



(1) 
(2) 



where px = \pt\, and the normalization is j d 2 pTfi(x,PT) = fi(%) and similarly for D\. Strictly speaking the Gauss 
widths (p^), {K\) could be flavor dependent, and x- or z-dependent, i.e., the Ansatze ([I] [2j in general do not imply 
a factorized x- or z- and transverse momentum dependence. 

Although in this way the unknown p^-depcndence is reduced to "one number," a Gauss width which "characterizes 
the pT-dependence" of a TMD, and although the Ansatz has not the correct large- pj asymptotics [9l|, it nevertheless 
proves to be useful in many processes sensitive to intrinsic transverse momenta j6ll - l63j . In fact, it was shown in [6l| 
that the Gauss model provides a useful approximation in many processes sensitive to intrinsic transverse momenta. 

The information on the Gauss model parameters in ([TJ [5]) used in many recent works was due to two independent 
studies of SIDIS data [H, [63[ . In the EMC data [32j on the azimuthal asymmetry Al^j in unpolarized SIDIS 
were used to fix the parameters in (JTJ [2]) - In (63j the widths were determined using (uncorrected for acceptance effects) 
mean values from HERMES on average transverse momenta {Phi) of produced hadrons. The results were found 



(Pt) 



0.25 GeV 2 
0.33 GeV 2 



m 
in 



62]. 

63], 



(K 2 T ) 



0.20 GeV 2 
0.16 GeV 2 



m 
in 



(3) 



which are consistent within the unestimated uncertainties involved in the analyses. Both studies were probably the 
best one could do at that time, but are not free of criticism. In [62] it was assumed the observable A^^f is due to the 
Cahn effect only, omitting effects of other TMDs. In [63| it was assumed the acceptance corrections were not large. 
In both works the very applicability of the Ansatz ([TJ [5]) was presumed. 

Meanwhile new data emerged, which improve the situation and allow to test the Gauss Ansatz in SIDIS more 
thoroughly than it was possible in [6ll463l ]. One of the aims of this study is to provide such tests. In Sect. [TT1 after 
convincing ourselves that the Gauss model passes the tests imposed by the new data, we will be in the position 
to update the information on the Gauss model parameters. We shall also discuss the EMC data on azimuthal 
asymmetries in the light of the updated information, and comment on the Cahn- and Boer-Mulders effect. 

Another purpose of this work is to test the Gauss Ansatz in the Drell-Yan process, and to determine the Gauss 
width of Also here we shall discuss azimuthal asymmetries in unpolarized DY, and comment on the Cahn- and 
Boer-Mulders effect, see Sec. Mil 

Finally, having established the applicability of the Gauss model in SIDIS and DY, we will address the question 
whether the descriptions of p^-dependences in the two processes are compatible. This is what one expects on the 
basis of the factorization approach and universality arguments, but a meaningful comparison of the effects requires to 
carefully take into consideration the different energies typically probed in SIDIS and DY, see Sec. lIVl Our conclusions 
are contained in Sec. [V] 
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II. INTRINSIC TRANSVERSE MOMENTA IN SIDIS 



HADRON PRODUCTION PLANE 



In this Section we study the distributions of transverse 
hadron momenta Phi. in SIDIS. Two effects play a role, 
namely intrinsic transverse parton momenta from the target 
which we model by (fT]), and transverse momenta the hadrons 
acquire in the fragmentation process which we model by ([2]). 
After a brief introduction to SIDIS, we shall study in detail 
data from CLAS and HERMES, and comment then on the 
Cahn effect at EMC and forthcoming experiments. 

Let P, I (V) and Ph denote (respectively) the momenta of 
the proton, incoming (outgoing) lepton and produced hadron. 
The relevant kinematic variables are q = l — V with Q 2 = —q 2 , 
W 2 = (P + q) 2 , x = Q 2 /(2P -q), y = (P ■ q)/(P ■ I), and 
z=(P-P h )/(P-q). 

In the following we are interested in the transverse momentum Pf,± of the produced hadron, which is defined with 
respect to the momentum of the virtual photon, see Fig. [TJ Assuming the Gauss model (JXJ [5]) the cross section 
differential in x, y, z, Pfij_ (but averaged over the azimuthal angle <p of the produced hadron) reads 




LEPTON SCATTERING PLANE 



FIG. 1: Kinematics of the SIDIS process Ip — > I'hX. 



d i cr uu (x,y 1 z 1 P hl _) _ An 2 a 2 s 



dx dy dz dP? 



h± 



J2e 2 a xf?(x)Dl(z)g(P hl _) 



(4) 



where the subscript "UU" indicates that both leptons and nucleons are unpolarized. The function G(Ph±) is given by 

1 



Q{P, 



exp 



ttk 2 ^(z) *~ \ K?r(z) 



p2 



k 2 t (z) = z z (p 2 T ) + {K 2 T ) , 



with the normalization J d 2 Ph±G(Ph±) = 1, while the cross section (j4]) is normalized such that 



d 3 auu(x,y,z) Aita 2 s 



dx dy dz 



Q 4 



l-y 



Y,e 2 a xfi{x)Dl{z) 



(5) 



(6) 



In the cross sections (j4} and (j6|) we neglect power suppressed terms of the order 0(M 2 /Q 2 ) where M is the nucleon 
mass. The neglected terms include purely kinematic factors, as well as a structure function (Fiju.l in the notation of 
[88j |) which has no partonic description. Depending on the kinematics these contributions of 0{M 2 /Q 2 ) need not to 
be small. We will recall this when necessary, but we will not need their explicit forms here. 

For later convenience we introduce also the notion of average transverse hadron momenta and their squares which 
are defined, and in the Gauss model given, as 



(Ph±(z)) 



Ph± dauujz.Ph^/dzdP 2 ^ 
dauu{z, Ph±)/dz dP^ 

P^dauu&P^/dzdP^ 
dauu(z,P h± )/dzdP^ x 



Gauss 



r y/ k t( z ) 3 



2 

4( z ) ) 



(7) 

(8) 



where (...) denotes average over Ph± ■ Analogously one could define mean transverse momenta as functions of other 
kinematic variables, for example (Phx(x)). Strictly speaking, in obtaining the Gauss model results in (JT1 |8j) we assume 
the Gauss widths to be flavor- and ^-independent. As we shall see below, these are reasonable approximations. 
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A. Lessons from CLAS and Hall-C 

In the CLAS experiment 33] the semi-inclusive 7r + electro-production off a proton target was studied with a 
5.75 GeV beam. Among others the following quantity was measured: 



R(P h ±) = 



_ d 4 a uu (x,y,z,P hl _)/dxdydzdPl^ 



d 4 auu(x, V, z, 0)/dx dy dz dP 2 _ 



= exp 



p2 



(9) 



at x = 0.24 and z = 0.30 for three different values of Q 2 . In the last step of ((9|) we assumed flavor-independent Gauss 
widths. For all three values of Q 2 the data are remarkably well described by the Gauss model [33j. In Fig. [5^ we 
show the data on the ratio © for the highest Q 2 — 2.37 GeV 2 , which is very well described with the parameter 



c t( z ) 



0.17 GeV 2 . 



(10) 



z=0.30 



The agreement of the data with the Gauss Ansatz is astonis hing . At lower Q 2 = 1.74 GeV 2 and 2 GeV 2 the situation 
is equally impressive, with somewhat lower values for Kj,(z) [331 ] . 

However, several reservations need to be made. First, in the CLAS kinematics [33[ the contributions of 0(M 2 /Q 2 ) 
we mentioned in the context of Eqs. ((4]) and ([6]) are not negligible, i.e., the description of CLAS data in Eqs. (l9l [T0|) 
effectively parametrizes also these contributions. Second, at z — 0.30 the measured hadrons are not only due to the 
fragmentation of the struck quark ( "current fragmentation" ) but can also originate from the hadronization of the 
target remnant ("target fragmentation"). The latter is described in terms of so-called fracture functions, which — 
while being a fascinating topic by themselves — are an undesired contamination from the point of view of DIS. 

Nevertheless, although one has to keep in mind these reservations, presently the 5.75 GeV beam CLAS data [33[ 
provide the best support for the applicability of the Gauss model in SIDIS. It would be desirable to solidify this 
observation with data taken at higher energies at CLAS12, HERMES and COMPASS. 

Another conclusion one can draw from the CLAS data [33|, modulo the above-mentioned reservations, is that 
the Gauss width (jxL) of fi(x,px) is only moderately ir-dependent, see Fig. which shows the average transverse 
momentum square (P 2 ^) of n + produced at z = 0.34 and Q 2 = 2.37 GeV 2 at CLAS [33[ as function of x. In fact, for 
0.2 < x < 0.5 we find (P 2 ±) = 0.15 GeV 2 within 20%, which is demonstrated by the shaded region in Fig. [2]d. 

In the Gauss model, the ^-dependence of the average hadron transverse momentum square (P 2 t (z)) in Eq. ([Sjl 
allows in principle to fix the Gauss widths (Kj) and (pj) of D\ and /f . Such data were presented in [33j, but here 
we will not use them for a quantitative determination of the Gauss model parameters. In fact, at the moderate beam 
energies for z < 0.4 (4lj the produced hadrons receive also contributions from target fragmentation, see above, and at 
still lower z "threshold effects" play a role [33]. Therefore we refrain here from using these data quantitatively, and 
will come back to them later for a qualitative comparison in Sec. Ill CI 
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FIG. 2: (a) The ratio R(Ph±) as defined in Eq. ([9| as function of the hadron transverse momentum square P 2 ±- The data are 
for 7r + from CLAS |33|]. The dotted line is an effective description in the Gauss model, see text, (b) The average transverse 
momentum square {P£±) of n + produced at z — 0.34 and Q 2 = 2.37 GeV 2 in SIDIS at CLAS [33| as function of x. The dotted 
line is an effective description in the Gauss model assuming the Gauss width of f" (a;, px) to be x-independent. This describes 
data within 20% in the region 0.2 < x < 0.5 as the shaded region shows. 
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Next we turn our attention to the Jefferson Lab data from the Hall-C collaboration 
were scattered off proton and deuterium targets in the kinematics 0.2 < x < 0.5, 2 GeV 2 < Q 2 < 4 GeV 2 



where 5.5 GeV electrons 
0.3 < z < 1, 



and tt with transverse momenta up to P^ ± < 0.2 GeV were measured. In spite of the narrow P^-range covered, 
the results on the differential cross-section (fi and E denote solid angle and energy of produced the hadron h with 



d 5 a(P, 



'hX Gauss 



d 5 a(0) e 



'hX 



dn dE dz dP 2 ± 



dfidEdzdP^ 



exp 



r>2 
_J_h±_ 

kL{z) 



(11) 



allow a valuable cross-check at (x) = 0.32 and (z) — 0.55 (t — p, d denotes the proton, deuteron target) 

As the kinematics is similar to CLAS, we expect the Hall-C data, which refer to (z) = 0.55, to be described in the 



Gauss model by the parameter K^{z) = 0.24 GeV . This value is taken from CLAS data on 



tf£i(*)> 



[33[ (see below, 



Fig. [5b). In this way we obtain a very good description of the P/^-dependence of the Hall-C data, see Fig. [31 

At this point it is interesting to stress that we assumed k\{z) = 0.24 GeV 2 to be the same for any hadron h from 
any target t. This means, we assumed flavor- independent Gauss widths (p 2 -,) and (-Kj>). Thus, although a dedicated 
analysis in 34| indicated a preference for slightly flavor-dependent Gauss widths, we conclude from Fig. [3] that the 
assumption of flavor-independent Gauss widths (v%) and (Kj) is reasonable — in the valence-x region for (z) = 0.55. 

To summarize, the data from Jefferson Lab [33,|34| are very well described assuming a Gauss distribution of intrinsic 
transverse parton momenta in the unpolarized distribution and fragmentation functions, and suggest a weak flavor- 
and x-dependence of the Gauss width (jpj) of fi(x,pr), and a weak flavor-dependence of the Gauss widths {Kj) of 
Df(z, Kt). However, one has to bear in mind that these data contain contributions from fracture functions and/or 
terms of 0(M 2 /Q 2 ), such that we shall limit ourselves here to these qualitative conclusions, and use for a quantitative 
analysis data from HERMES [H, see Sees. QTB] and QTC] 
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FIG. 3: The differential cross-section d 5 <r/df2 dEdz dP%± for pion production off proton and deuterium targets at (a;) = 0.32 
and (z) — 0.55 as function of P^± from Hall-C [3J|. The theoretical curves are from the Gauss model with the Gauss width 
fixed from CLAS [33l ]. The overall normalization of the cross sections is fixed by hand. 
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B. Testing the Gauss model at HERMES 

Assuming the Gauss model (Q] [2]) and flavor- and x-independent Gauss 
widths, we obtained for the mean transverse hadron momenta and their 
squares the results quoted in {[71 [HJ) - In particular, these two quantities are 
related in the Gauss model according to 

(P^(z)) 2 = \{Plx_{z)). (12) 

Meanwhile data from HERMES on SIDIS of electrons or positrons with 
-Ebeam = 27.6 GeV off a deuteron target allow to test the Gauss model 
prediction (JI2J) in the kinematics Q 2 > lGeV 2 , W 2 > 10 GeV^w < 0.85, 
z > 0.2 and 0.023 < x < 0.4 with (Q 2 ) = 2.4 GeV 2 , (x) = 0.09 jsEHl. 

In Fig. @] we show (Ph±{z)) as function of z from [3!| (triangles). These 
mean values refer to pions and kaons and were not corrected for acceptance 
effects. We compare them with \\^{Ph±_{ z )) 1 ^ 2 using the data from [35[ 
on (P 2 ± (z)} for positive pions. At HERMES the SIDIS events are subject 
to the cuts 0.2 < z < 0.7, as indicated in Fig. 01 and in this range 7r ± and 
K + have very similar {P% ± (z)) H|. We shall come back to this observation 
in the next section. Let us add that the actual aim of [35| was to study 
nuclear Ph± -broadening effects, the so-called Cronin-effect, and various 
nuclear targets were used besides deuteron, namely He, Ne, Kr, Xe. 

We conclude from the exercise in Fig. [4] that the SIDIS data from HERMES support very well the Gauss model 
relation (|12p. even though one should keep in mind possible effects due to acceptance corrections in the case of the 
(Ph±(z)) results fro m 39l| . Of course, data on the cross section differential in P^± from HERMES, of the type as 
presented by CLAS [33j , are required to provide fully conclusive support for the Gauss model. However, we still may 
view the picture in Fig. [4] as an encouraging indication in favor of the applicability of the Gauss model in SIDIS at 
HERMES. 



C. Determining parameters from HERMES, and cross check with CLAS 

The HERMES data [35J allow two further insig hts. The first is that the (P* ± (z)) for 7r + , 7r , K + are very similar 
in the SIDIS region, i.e., below the cut z = 0.7 which excludes exclusive effects at HERMES, see Fig. [SJi. Thus, there 
is no evidence of a strong flavor-dependence of the Gauss widths (p\) and (Kj,) at HERMES. The second insight is 
that the data in Fig.^ allow to fix the Gauss widths (pi^) and (.Kj.). A fit in the region 0.2 < z < 0.7 yields 

(p 2 T ) = (0.38 ± 0.06) GeV 2 , 

(Kl) = (0.16 ±0.01) GeV 2 , (13) 

with a x 2 P er degree of freedom of 0.44. The best fit and its 1-a region are shown in Fig. [SJi as (respectively) dotted 
line and shaded region. 

The values of the new fit in Eq. (fT3|) are in good agreement with the results from [H,|63| quoted in Eq. ([3]), especially 
if one recalls that those numbers have unestimated systematic uncertainties. This is the main point, in which our 
new result (fT3| constitutes an improvement over the old numbers ([3]): it has a well-estimated uncertainty. 

In any case, the good news is that the numbers for the Gauss widths from the works (621 . [63| are confirmed 
numerically. This means that, in this respect, phenomenological results for fixed target experiments obtained on the 
basis of those numbers remain valid. 

Finally, let us turn back to the CLAS data [33[ on the z-dependence of (P%^(z)) which refer to (Q 2 ) = 2.37 GeV 2 
which is similar to HERMES, but (x) = 0.27 which is substantially higher than HERMES. In Sec. lHAl we refrained 
from using these data to determine quantitatively Gauss model parameters for reasons discussed there in detail. At 
this point it is instructive to compare the fit result obtained from HERMES to the CLAS data, see Fig. [S}d. Clearly, 
in the region above z > 0.4 where the produced pions are predominantly from current fragmentation, we observe a 
good agreement. This indicates that it is the same non-perturbative mechanism which generates intrinsic transverse 
momenta in the two experiments. 
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FIG. 4: Transverse momenta of hadrons 
in SIDIS off deuterium at HERMES vs. z. 
We compare (Ph±(z)) (triangles) from 3^ ]. 
with | {Pi ± (z)) 1/2 (squares) from [H]. 
In the indicated SIDIS range of HERMES 
these quantities are predicted to coincide 
in the Gauss model, see Eq. (I12[) . 
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FIG. 5: (a) (P% ± (z)) of hadr ons in SIDIS off deuterium at HERMES vs. z from [3a|. The dotted line and the shaded region 
are the best fit and its 1— <x region from Eq. JT3J, see text, (b) (P k 2 j_(?)) of ir + in SIDIS off proton at CLAS vs. z [33t ] ■ Dotted 
line (shaded region) are the best fit (its 1— a region) to HERMES data from Fig. [5^,, see text. 



D. Cahn effect at EMC 



In SIDIS of unpolarized leptons off unpolarized nucleons the cross section differential in the azimuthal angle of the 
produced hadrons, see Fig. Q] is — for purely electromagnetic interactions — given by [8^ | 



d 5 a{x,y, z,P) 



h±, 



dx dy dz dP% ± <. 



d i a uu (x,y,z,P hl _) 
(2tt) da; dy dz dPj 



(2 - 

l + (l-y) 2 



cos(</>) A c °Ij 4 



y 



cos 20 



l + Cl-y) 1 



cos(20) A™ v 



(14) 



where d 4 arj(j(x, y, z, Ph±) denotes the differential cross section d 5 auu(x, y, z, Ph±, 4>) after averaging over <f), as defined 
in Eq. (TJ}. A^f/ and A™^ 20 in (fH]) depend on x, z, Phi. and are ratios of adequately defined structure functions. 



The meaning of the subscript "UU" is as in the superscript recalls the type of 0- modulation. 

At low Ph±, the observable A^^f is suppressed by 1/Q, and in this sense a "twist-3" effect. At present it is 
not clear whether there is factorization in SIDIS at subleading twist [9(J HH . This cos 0- modulation is sometimes 
referred to as the Cahn effect. In Ref. |9f it was shown that the existence of intrinsic transverse parton momenta 
in the unpolarized distribution and fragmentation functions can generate such a modulation. Later it became clear 
that, if one assumes factorization and restricts oneself to a lowest order QCD ( "tree- level" ) description, there are four 
different contributions to this asymmetry from several TMDs and -ftf^-dependent fragmentation functions [T7l - fl9l l88j . 
More precisely, in the notation of [H| the asymmetry is given by 



a cos (J) _ r UJJ 

A uu - Ip 



cos <f> 



UU,T 



Ft 



UU 



2M 



Fuu' 



K T ■ P h ± 
m h Ph± 



C 



fiDi 



xhH^ 



m h D x 
~M fl — 



PT ■ PhA 

MP hX 



m h h ±H\ 



(15) 



where denotes the hadron mass, and the meaning of the convolution symbol is, for generic functions / and D, 



• fD 



J d 2 p T J d 2 K T 6 {2) (zp T + 



K T - P hl _) 



Y,^xf a (x, PT )D a (z,K 7 



(16) 



The original "Cahn-effect-only" description of ^4^^ 0] can be "rederived" from the TMD formalism [!' 



11 



under two conditions. The first is that quark-gluon-quark correlators, which can be separated off by exploring QCD 
equations of motion, can be neglected with respect to quark-quark correlators. The quality of such approximations was 
discussed in [92H95l |. One of these approximations is xf ±q ~ ff, and is supported by results from the bag model [9fj| . 
The second condition is that in the experiment the type of hadron h is not detected, i.e., the observed hadron with 
momentum specified by z, 4>, Ph± is a sum over pions, kaons, etc., and in such sums the Collins effect tends to cancel. 
Strictly speaking, an exact cancellation occurs only if one weighs the Collins function H x adequately, integrates over 
< z < 1, and sums over all hadrons (97[. However, string fragmentation models [H| and phenomenology [99Hl02j ] 
indicate that cancellations occur in practice already if one sums over charged hadrons only, which are mostly in 
experiments such as EMC, where the asymmetry A^^ in charged hadron production was measured (3lL l32|. 
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In fact, the EMC data on A 



COS 

UU 



13 ll. 13211 were used 



in [62| to determine in the "Cahn-effect-only" approxi- 
mation the Gauss model parameters quoted in Eq. <j3j) . 
It is therefore instructive to check, whether the revised 
numbers from (|13|) still give a good description of these 
data. 

In the "Cahn-effect-only" approximation, i.e., ne- 
glecting in (fl"5"j) the "pure twist-3" tilde-functions and 
the Collins effect, we obtain for the asymmetry with 
the Gauss model JTJ [5]) the result 



cos((j))(z) 



A 



U U. Cah 



/7F (p; 



(Q)y/z*(j? T ) + (K*) 



(17) 



where we assumed flavor-independent Gauss widths. 
In the EMC 31| experiment 280 GeV muons were 



scattered off protons and the covered kinematics was: 
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-0.1 
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EMC data $ 
Cahn effect 



> 0.15, Q 2 > 



10 GeV 2 with about (Q) 
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FIG. 6: Azimuthal asymmetry in charged hadron pro- 

duction vs. z. The data are from the EMC experiment [311 ]. 
The theoretical curve is the "Cahn-effect-only" approximation 
for this observable, which is justified under certain assumptions 
(see text), using the Gauss model with parameters fixed from 
HERMES, Eq. Q3). 



V < 0.8, 

4.8 GeV, and 160 < W 2 /GeV* < 360. Also the cut 
Phi_ > 200 MeV was imposed which is ignored in (|T7|) 
for simplicity. Numerically it has a negligible effect. 
In this kinematics we obtain with the Gauss model pa- 
rameters inferred from HERMES, Eq. (fTBl . the result 
shown in Fig. |6l We observe a very good agreement. 

At this point, this not only demonstrates the com- 
patibility of the EMC and HERMES data. Since we 
have fixed the details of the Gauss model in an inde- 
pendent experiment, the excellent agreement we ob- 
serve in Fig. [5] can also be read in opposite direction. 

The approximations we used in order to relate A™^ 1 to the Cahn effect, namely the cancellation of the Collins effect 
in charged hadron production and the neglect of tilde-terms, are justified — within the experimental error bars and 
the theoretical uncertainty of our study. In this sense, the EMC data on A^^ support the "Wandzura-Wilczek-type" 
approximations discussed critically in [92T - |95| . (However, we shall come back to this point at the end of Sec. IIVI ) 

To draw an intermediate summary, in the Sections III AMI Dl we have seen that data from EMC [HI HH , Jefferson 
Lab [111, [34[ and HERMES 35] support the Gauss model with flavor- and x- or z-independent Gauss widths. In the 
future one may need to refine the description by allowing for flavor- and x- or z-dependent Gauss widths, when more 
precise data will make it necessary to introduce and possible to constrain further parameters. 

The Gauss model has an important principle limitation though. It may work only if the transverse momenta of 
the produced hadrons are of the order of magnitude of the hadronic scale, i.e., much smaller than the hard scale in 
the process. This condition is fulfilled in the case of EMC, Jefferson Lab, and HERMES data [3ll - l35| where typically 
(Ph±) — 0.4 GeV <C (Q) = 2-5 GeV in these experiments. When transverse hadron momenta become substantially 
larger than the hadronic scale or even become so large that they set the hard scale in the process, one does not 
deal with non-perturbative intrinsic pt anymore, but can apply perturbative QCD [8]. SIDIS data from high energy 
experiments, for example E665 at Fermilab [3(| or ZEUS at DESY [53, HH, are sensitive to such perturbative pr- 
effects. In practice, in order to describe these data it is necessary to in clude both, perturbative and non-perturbative 
effects. We refer to the pioneering study [l(|, see also the recent work llO 
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E. Cahn and Boer-Mulders effect, and new data 



It would be interesting to repeat the analysis presented in the previous section with recent results on A^y from 
Jefferson Lab |H, H3] (final), HERMES and COMPASS jM SI (preliminary data). Since (Q) in these experiments 
is smaller compared to EMC, one would expect this subleading- twist observable to be there larger than at EMC. 

A careful comparison of the Cahn effect prediction for A™^ to these new data will be instructive for two reasons. 
First, at Jefferson Lab, COMPASS, HERMES the hadrons arc identified, i.e., the Collins effect in (fT5j) does not cancel 
out — in contrast to EMC, where we could argue it may (at least approximately) cancel out. Second, the new or 
forthcoming data are far more precise compared to EMC, such that deviations of data from the "Cahn-effect-only" 
approximation will have a chance to pop up more easily, and provide insights on quark-gluon correlations. In view of 
the preliminary status of the HERMES and COMPASS data [48|, |49(, however, we shall not pursue this study here. 

Let us also comment on the Boer-Mulders effect responsible for the cos(2(/>)-modulation in (fH| and given by 



cos 20 



A cos 20 _ r UU 

uu ~ Fan: 



O 



(Ml 



pcos 20 



2(K T ■ P h ±)(pr ■ Ph±) - (K T 



Mm h P 2 hl 



H 



(18) 



is one of the 



with Fuu,t and the convolution integral as defined in (jT5l fT6]) . The Boer-Mulders function \li 
so-called "naively time-reversal odd" (T-odd) TMDs (23U25] . The observable was discussed in |64471j . 

In Eq. (1181) we have indicated the power corrections modulo which the factorization theorem is formulated, e.g. 
pi. [Til Il2l l26j| . Only if such corrections are sufficiently small the factorization approach is justified, and can develop its 
predictive power with universal non-perturbative objects [27] . though the concept of universality had to be extented 
in order to accommodate T-odd TMDs [23U25I ] . which we will discuss in more detail below in Sec. MI Cl 

In general such non-factorizable corrections are theoretically not under control and must be excluded experimentally 
by studying the Q-behavior of the cross section or asymmetry. In the case of A™^ 2 ^, however, the parton model 
provides a way to estimate one of the possible power corrections in ([TBI) . Namely, the Cahn effect [9j yields 



*cos 20 



(P 



(Q 2 ) z 2 (p 2 T ) + (K T ) 



(19) 



if we assume the Gauss model with flavor-independent Gauss widths. We stress that this is only one among many non- 
factorizable power-corrections to the cos(2</>)-modulation. Ideally, one should choose the cuts in Q in the experiment 
such that one does not feel the Cahn- (fT!)|) or any other 1 / Q 2 -power correction in (fT5|) . 

In practice, this is difficult. At HERMES and COMPASS the DIS cut Q 2 > 1 GeV 2 is imposed which leads to 
(Q 2 ) ~ 2.5 GeV 2 . At Jefferson Lab similar (Q 2 ) (but lower W) is achieved, in spite of lower beam energies, because 
there one can afford to choose the cut Q 2 > 2 GeV 2 [HJ thanks to the high luminosity. In such kinematics, at let us 
say z = 0.5, the Cahn-power correction (TT9l is, with Gauss widths from (fT5|) . about 



a cos 20 
A UU,Cahn 



-0(5%) 



(20) 



Let us compare this number to what one may expect from the leading twist Boer-Mulders contribution in (I18[) . If we 
generously allow the Boer-Mulders function to saturate its positivity bound ^-\h^ a (x,pT)\ < fi(x,pr) [21], then 



^cos 20 
UU 



(x,z)\<C 



WBM ji \J^i I 

m h 



_L(l/2)o 



(21) 



where cjbm = %(Kt ■ Ph±)(j>T ' Ph±) ~ (Kt 'Pt), with the "hat" denoting unit vectors and the transverse (l/2)-moment 
defined as H^ 1 ' 2 ' a (z) = J d 2 KT^H^- a (z, Kt)- The numerical factor ctQ — 0(1) contains the dependence on the 

model for the transverse parton momenta. Phenomenological studies indicate that |i? 1 L( ' 1 ^ 2 '' a /_D°| = 0(10%) around 
z = 0.5. Thus, the leading-twist Boer-Mulders contribution in (TT51) can be expected to be at most 

\A C ut4 M \<G(10%). (22) 

Hence, the power correction estimated on the basis of the Cahn effect in ([20]) is substantial in this kinematics. Studies 
of the Boer-Mulders effect in SIDIS [63 - [70l | and first analyses of preliminary data [7l| yield similar results. We 
conclude therefore that, unless it will be possible to suppress such power corrections by applying sufficiently large cuts 
in Q, or separate experimentally the Boer-Mulders part and power corrections by their different Q-behavior, the task 
of gaining insights on the Boer-Mulders function from present SIDIS data will be highly demanding — and feasible 
only on a basis of a thorough quantitative understanding of intrinsic transverse momentum effects in unpolarized 
distribution and fragmentation functions, which is what this work aims at. Interestingly, as we shall see in Section [TTT1 
the situation is less demanding for the available DY data. 
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III. INTRINSIC TRANSVERSE MOMENTA IN DRELL-YAN 



In this Section we discuss the DY process |104| . i.e., the inclusive production of large- invariant-mass /i + /i pairs 
in hadron -hadron collisions hih 2 — > fJ, + fi~X, first observed in 1970 at the Brookhaven AGS |105j . for reviews see 
- IllOj . Let pi^ and k\ t 2 denote (respectively) the momenta of the incoming hadrons hi 2 and the outgoing lepton 
pair. The kinematics of the process is described by the center of mass energy square s, invariant mass of the lepton 
pair Q, rapidity y or the Feynman variable xp, and the variable r which are defined as 

s = (Pi+P2) , Q ={ki+k 2 ) , y= -In T — —— = -ln— , 

2 Pi ■ («i + k 2 ) 2 x 2 

Q 2 

Xp = X\ — X2 , T = = XiX 2 . (23) 

s 

The three-momentum of the virtual photon q = ki + k 2 can be decomposed in the center of mass frame of the incoming 
hadrons into a longitudinal and transverse component with respect to the collision axis as q = {(IltQt)- I n this frame 
y = \ ln(Q + Pl)/{Q - Pl) and x F = 2q L /^/s. 

In the parton model Xi denotes the fraction of the hadron momentum pi carried by (respectively) the annihilating 
parton or anti-parton, and the cross section differential in Xi, X2 and qT = \qr\ is given by (we write fi(xi) = f^/ hi (xi) 
for brevity, i.e., the parton distributions in the possibly distinct hadrons hi are unambiguously labelled by the Xi) 

d 3 auu kira z 



e 2 a J r d<j>J d 2 PlT J d 2 p 2T 5 {2 \p 1T +p 2T ~ q T )fi{xi,p 1T )fKx2,P2T) ■ (24) 



dxi dx2 qr dqr 9Q 2 

a—u,d,u,... q 

For an accurate, quantitative description of the ^-dependence of the DY cross section it is necessary to go beyond 
the parton model expression (|24p and use the Collins-Soper-Sterman formalism [12j . 

However, it turns out that a very useful effective description can be obtained by assuming the Gauss model ([TJ. 
Inserting (p} into (f24| . and assuming that the Gauss widths {{p\t) in hadron 1, and (p 2T ) in hadron 2) are flavor- 
independent, yields 

a — a — a — = ~^Fa 1^ e afi{xi)h{X2)2q T 2 , K DY = (p 1T ) + (p 2T ) . (25 

dxi dx 2 dqr 9Q , _ K n v 

a—u^d,u,... ^ 1 

We define the mean transverse momentum (qx) of the muon pair, and the mean transverse momentum square (q^) as 

2 / q^dauu/dqrX Gauss 2 

Wr> = < -h -n > = k dy , 

\ dauu/dqr / 

/ qT dauu/dqr \ Gauss \pK !~2 ( . 

{QT) = { dauulMr / = -V^" < 26 > 
When integrating ()25[) over qx the transverse momentum model dependence drops out and one obtains 

d 2 auu ^iru 2 



E 4/i°(zi)/?(*2) ■ (27) 



dx\ dx 2 9Q 2 

a—u,d,u,... 

Using for /f (x) leading order parametrizations from DIS, the leading order QCD formula (j2"T|) underestimates data 
by a ("K"-)factor ~ (1.1-1.7) depending on Q 2 (and moderately also on 2:1,2) [l06t | . 

One gets the overall normalization of the DY cross section correct only by including higher order QCD correc- 
tions. Throughout we shall work here in leading order. It is a remarkable observation that, whenever corrections 
to leading order DY cross sections were considered, absolute cross sections were found to be considerably altered 
(typically enhanced), but ratios of observables in DY were found to be only moderately altered. This is in par- 
ticular the case for collinear double spi n asymm etries, for example the double transverse spin asymmetry Att oc 



ymme 

huh. 
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A. 



Testing the Gauss model in DY 



The Gauss model unambiguously connects the average lepton pair transverse momenta and their squares. In fact, 
from (f2"rJ)) we obtain 



which is a parameter-free prediction of the Gauss model, c.f. Eq. (fl"2)) in SIDIS. Data allowing to check (|28|) were 
reported in [52j, and are shown in Fig. [7j Several interesting conclusions can be drawn from Fig. [7j First, it 
demonstrates that in the range 0.2 < xp < 1 covered in the experiment [HJ the relation (|28p is satisfied within the 
accuracy of the data. Second, it demonstrates that for 0.2 < xp < 0.9 both (g^) and (gr) are xp -independent within 
error bars. Therefore, we conclude that for xp < 0.9 the Gauss model with ^-independent Gauss widths provides an 
excellent description for the ^-dependence of the DY process, which is useful at least for the description of data with 
an accuracy comparable to that of [52| . More precise data may reveal deviations from the Gauss Ansatz. 

It should be n oticed that the parton model is not adequate at large xf, where higher twist effects are expected to 
be relevant [ll7j . Interestingly, the Gauss model relation (|2"5|) holds even in the region xp > 0.9 within error bars, 
see Fig. [7J 

In Ref. [l!8| (gf.) was measured for DY lepton pairs induced by various hadrons on a platinum target. Within 
experimental accuracy it was found that (q T )n- Pt ~ (^T^+ft j e ^ c - Since these reactions arc dominated by differently 
flavored valence quarks, one may conclude that it is indeed approximately justified to us e flavor-independent Gauss 
widths. It is also important to notice, that the dependence on the nuclear target [119] is moderate for transverse 
dilepton momenta qt < 3 GeV, i.e., nuclear effects are not relevant. 

The comparison in Fig.[JJis very instructive, and allowed us to draw many interesting conclusions. Yet, it is not a 
proof that the Gauss model really works in DY. For that it is necessary to consider data on cross sections, which we 
will do in the next section. 



(Qt) 



2 \ Gauss 4 




(28) 




t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 
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FIG. 7: The mean dimuon transverse momentum square (g^) vs. if as measured in the Fermilab E615 experiment [Kg . The 
data points for {q^} are marked by circles, the data points for ^{qr} 2 are marked by triangles. Both quantities are predicted 
to be equal in the Gauss model, see Eq. (|28[). which is the case within the statistical accuracy of the data. 
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B. QT-dependence of the cross section 



The qx- dependence of the cross sections allows to fix the parameters of the Gauss model. It is convenient to rewrite 
the cross section (f2"5"j) as (here ajjjj denotes the total (^-integrated unpolarized DY cross section) 

dauu(q T ) 2q T exp(-q T /K% Y ) „dcr l/[ /(g T ) 2 Aa vu exp(-q^/ k 2 dy ) 

j = a uu 2 ' or y m = ~t~a 2 ■ 

aqr k dy a A q y s clxf ^ dy 

Fig. M (left) shows FNAL-288 data on the qr-dependence of the invariant differential cross section Q^f- for 
production from scattering a 300 GeV proton beam off platinum at (y) — 0.21 for two different bins in Q j50[. Clearly, 
the Gauss model provides a good description of the data for 

4r h =2<^ r ) = 1.4GeV 2 , (30) 

where ip 2 N j) denotes the mean parton momentum square in the nucleon. 

The limitations of the Gauss Ansatz become more evident by using a logarithmic scale for the cross section. We 
do this in Fig. [5}d that shows FNAL-E615 data [55| for the differential cross section for (i + fi~ production from 
scattering a 256 GeV tt~ beam on a tungsten target for < xp < 1 and 4.05 GeV < Q < 8.55 GeV. The Gauss model 
is applicable for q? < 3 GeV with the parameter 

4r h = (p^) + <P 2 wT ) = 1.7GeV 2 ) (31) 

where (p 2 y) denotes the mean parton momentum square in the pion. 

The result (foTj) corresponds precisely to the experimental average (q 2 -) = (1.71±0.02) GeV 2 [H)]. This is so because 
the differential cross section in the region < qr < 3 GeV (where the Gauss Ansatz is applicable) is about two orders 
of magnitude larger than in the region qx > 3 GeV where it is not. Therefore the mean transverse momentum square 
is dominated by that region of qx where the Gauss model works. This means that in practice, if the Ansatz works, 
the respective Gauss widths can be read off directly from the experimental results for (q^). 

Since the kinematics is comparable in the FNAL-288 and FNAL-E615 experiments, we can use ([3U1 l3"Tj) to deduce 

(pl t ) = 0.7 GeV 2 , (pl T ) = 1.0 GeV 2 at ^ - 23 GeV. (32) 

Thus, transverse parton momenta appear to be larger in the pion than in the nucleon, as already noticed in (53j . The 
nucleon width in (1321) is in good agreement with the value (p 2 NT ) = l//? 2 = 0.64 GeV 2 obtained at i/s ~ 20 GeV in 
611. 




FIG. 8: Left: The invariant differential cross section for p Pt — > X at (y) — 0.21 for two different Q-bins from 

FNAL-288 [Ij. The Gauss model, Eq. (J23), provides a good description of the data for (gf.) = 1.4 GeV 2 . 

Right: The differential cross section for it~W — > y, + fi~X from FNAL-E615 Here the Gauss model, Eq. (f29)) . provides 
a good description of the data up to qr < 3 GeV with {q^} = 1.7 GeV 2 . 
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Cahn and Boer-Mulders effect in DY 



Next we turn our attention to the azimuthal dependence of the DY cross section in unpolarized hadron collisions, 
which was measured in var ious experiments [5ll457j | . The general expression for the angular differential cross section 
in the Collins-Soper frame [120l ] is traditionally written as 



a dfi 



1 



1 + A cos 2 9 + [i sin(26») cos < 



— sin 
2 



cos (2(f)) 



(33) 



4tt A + 3 

In collinear QCD factorization to 0(a s ) the coefficients A, fi, v are expected to obey the Lam- Tung relation [l2l[|l22| 



2u + A = 1 , 



(34) 



which is preserved as a good approximation at 0(a 2 s ) [l23| . But in experiments with pion beams on nuclear targets 
the relation was found to be strongly violated: while A = 0(1) was observed, large values v ^ were found 
[5ll - [55j . (However, in proton-proton and deuterium-proton collisions [H, H3| the relation (1341 was found satisfied.) 

An attractive explanation for the violation of the Lam- Tung relation ([34| is provided in the framework of TMDs 
f2~oj in terms of the Boer-Mulders function hj~ [lj|. In this approach v is related to ^ a e^h^ q (x\)h^ (X2) 20]. 

It is of interest to learn about the Boer-Mulders function from DY and SIDIS because this T-odd TMD was 
predicted to have unusual universality properties. On the basis of time-reversal arguments it was predicted [24j ] that 
h± in semi-inclusive deeply inelastic scattering (SIDIS) and in the Drell-Yan process (DY) have opposite signs, 



hi(x,p T ) D is = -hi(x,p T ) DY 



(35) 



Analogous relations are expected to hold also for the Sivers function and other T-odd TMDs [2J,l25|]. The experimental 
check of such universality relations for T-odd TMDs would provide a thorough test of our understandin g of the 
factorization approach to transverse momentum dependent processes in terms of pr-dependent correlators [24l - l27| . 
Perspectives to test the universality relation for the Sivers function where discussed in [72Tj80| . and for the Boer- 
Mulders function in |77H83j . 

As we have seen in Sec. Ill El the extraction of the Bocr-Muldcrs function from SIDIS is hampered by substantial 
power-corrections due to the Cahn-effect. It is interesting to ask whether the same difficulties occur also in DY. In 
order to address this question a quantitative understanding of intrinsic px-eSects in DY is necessary, and on the basis 
of the results from Section lill Bl we are prepared to have a look at th at. 

As in SIDIS, in the LO QCD ("tree level") approach at low <jt [llOl ] the coefficient of the cos</> modulation in (j3"3"| is 
suppressed by one power of the large scale Q, but v which is related to the Boer- Mulders effect is leading twist. And, 
as in SIDIS, the Cahn effect generates a l/Q-contribution to /i and a l/Q 2 -power-correction to v. These contributions 
are given by 



A*Calm 



^Cahn 



= A 



qr_ 

Q 



(P?T> 



(Par) 



Wit) + (pIt) 



A 



q 2 r 



(Pit) 



(Par) 



(P?T> 



+ (pIt) 



where 



.4 



1 + 



4(p 



It) (P2T 



2Q 2 



(P?t) + (P2T 




(36) 
(37) 

(38) 



Notice that in both cases the effect vanishes, if the Gauss widths are equal. Which does not mean that if vanishes, 
if the colliding hadrons are the same. For example, in proton-proton collisions the Gauss widths of quarks and 
antiquarks enter, which could be different. It is true that in Sec. [TT] we saw no evidence for a flavor dependence of 
the Gauss widths, but we should keep in mind that in SIDIS at Jefferson Lab and HERMES sea quarks do not play 
a dominant role. 
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FIG. 9: The coefficient fi and v in the azimuthal distribution (|33[1 of the dileptons in DY as function of the dilepton transverse 
momentum qr in the Collins-Soper frame. Data are from the Fermilab E615 experiment [551 ] . The curves represent an estimate 
of the Cahn effect contribution to these coefficients, see text. 



With the results on the Gauss widths for pion and nucleon inferred in Eq. (]32l) we obtain the estimates for the 
contributions of the Cahn effect to the coefficients fi, v shown in Fig. in comparison to Fermilab E615 data taken 
from 7r~-nucleus collisions [55j . According to the convention what is hadron 1 and hadron 2 in (55j we have to identify 
(p\ T ) = (pl T ) and (p 2 2T ) = (p 2 NT ) in Eq. ([32]). 

We see in Fig. |9K that the Cahn effect is able to account partly for the power-suppressed coefficient \i of the 
cos (^-modulation in (f3"3"]). At this point we should not be worried too much about the fact, that in SIDIS at EMC 
the Cahn effect was able to account for the entire cos (^-modulation. Here we deal with a different kinematics and 
TMDs replacing the role of fragmentation functions. (For the latter reason one also cannot conclude anything on the 
WW-type-approximations, see Sec. Ill DP Also we have to keep in mind that the result for ^ is strongly sensitive 
to the difference of the Gauss widths for pion and nucleon, which we estimated crudely in Eq. ([32]) . In fact, if 
one was interested in that, one could describe the /^-coefficient entirely in terms of the Cahn effect with the choice 
(p 2 T ) = 1.3 GeV 2 and (p 2 NT ) = 0.3 GeV 2 , but we shall refrain from doing this here. 

What is important at this point is the observation, that the Cahn-effect-power-correction to the ^-coefficient is 
negligible. We conclude from this exercise that the present DY data provide a "safer" way of accessing information 
on the Boer-Mulders effect, in the sense that they are far less sensitive to power-corrections as compared to available 
SIDIS data. In practice, of course, data from both processes need to be explored in order to test the univ ersal ity 
prediction (JSSJ). Studies of the Boer- Mulders effect were presented in SIDIS and DY [zHHl, see also [HI]. 
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IV. ENERGY DEPENDENCE OF INTRINSIC TRANSVERSE MOMENTA IN DY AND SIDIS 



In Eqs. (|30[ I31|) we were able to combine information on the ^-dependence of DY cross sections from different 
experiments to arrive at (|32l) . because of comparable kinematics. In general experiments performed at different 
energies have to be compared with care. In DY the mean lepton momentum square {q 2 -.) is energy (s-)dependent. 
This depend ence is different for i = ttN or pN induced Drell-Yan. For 50 GeV 2 < s < 600 GeV 2 it can roughly be 



described as 



(Qr( s ))i ^A. + B.s 



A^n = (0.59 ± 0.05) GeV 2 , B„ N = (2.8 ± 0.2) • 10" 3 , 
A pN = (0.52 ± 0.11) GeV 2 , B pN = (1.4 ± 0.2) • 10~ 3 . 



(39) 



This is just one way of parametrizing the s-dependence. 
The dat a ar e compatible also with a linear in \fs increase 
of (qj,) [125| . while in QCD one would rather expect a log- 
arithmic increase. But in a limited s-range an effective 
parametrization of the type (|39[) works reasonably well. In 
any case, (qf.) increases with energy, which reflects the trans- 
verse momentum broadening due to gluon radiation [l2j |. 

Fig.[TUlshows (qj,) in irN induced reactions and the original 



fit from [125[. The values of (q^nN deduced from Figs.[7]and 
[8] are shown for comparison but not included in the fits. As- 
suming the Gauss model (|26|) we obtain from ([39]) the follow- 
ing effective s-dependence of the intrinsic transverse parton 
momenta in the hadron h: 



(Pt(0)) 



C h 



04(0)) = 0.3 GeV 
C h = 10" 3 



2.1 for h = 7T, 
0.7 for h = p. 



(40) 
(41) 
(42) 



At this point two questions arise. First, are transverse 
parton momenta in DY and SIDIS compatible? Second, if so, 
is there any indication of transverse momentum broadening 
in SIDIS too? 

Concerning the first question, it is encouraging to observe 
that flU]) "predicts" for HERMES (s = 52 GeV 2 ) the result 
(Pt ( s ) ) | hermes = 0.34 GeV 2 which is within the error bars 
of the Gauss width in (fTBf . Probably it would be more consis- 
tent to use in SIDIS the photon-hadron center of mass energy 
square W 2 instead of s. However, in view of the uncertainties 
in (Tf"3"|) and (f4"01 — I4"2l this is numerically of little relevance. 

Concerning the second question, let us compare the mean 
square transverse momenta (P 2 ±(z)) in SIDIS from various 
experiments at a common value of z, let us say 0.5 < z < 0.6 
(as different ranges are covered we cannot compare averages 
over z). For Jefferson Lab we use 

( P h±( z )) = 0.24 GeV 2 at 
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FIG. 10: Mean dimuon transverse momentum square 
(qx) as function of the center of mass ener gy square, s, 
in ir~ N induced Drell-Yan. Following |l25l |. 
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FIG. 11: Mean square transverse momenta (Ph±{ z )} in 
SIDIS around 2 ~ 0.5 as function of s from Jefferson Lab 
UlSl. HERMES [U, COMPASS HJ. 



z = 0.55 from CLAS [33|, which describes well the Hall C data 
[13, see Figj3 For HERMES we take (P^A Z )) = 0.27 GeV 2 
at z = 0.52 [H. For COMPASS we use (P h ±{z)) = 0.55 GeV 
at z — 0.566 from [40J which we convert by means of (fL2|) . 

Fig. 1111 must be interpreted with care. The shown (P 2 ±(z)} were obtained in different ways, span a small s-range, 
and have systematic uncertainties except for the HERMES value. For a conclusive comparison acceptance corrected 
data are needed from all experiments. Nevertheless, we see a tendency for an increase in s with a slope which is 20 % 
lower than (|4"2]) . see Fig. [TTJ We made no effort to estimate the uncertainty of (|42|) but it is presumably not smaller 
than 20%. So the s-slopes in SIDIS and DY are compatible. 
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Notice that in SIDIS also the "broadening" of the width of the transverse momenta in the fragmentation function 
contributes, and we above tacitly assumed that both (p|n) in fi(x,pr) and (Jf|) in Df(z, Kt) increase with s at a 
similar rate. In fact, would we extrapolate from HERMES at s = 52GeV 2 , Eq. {T3J), assuming that {p\) increases 
with s at the rate (g2J but keeping (K%) at its initial value in (fI5J). we would "predict" (P, 2 ± (z)) ~ (0.30 ±0.03) GeV 2 
&tz = 0.566 for COMPASS which strongly underestimates the measured value (P^z)) ~ 0.39 GeV 2 0. Thus, the 
result in Fig. [TT] also indicates {Kj) -broadening in the fragmentation function. 

Finally, let us comment on EMC data on A^^. In Sec. HID I we have seen that the Cahn effect can explain these 
data with the Gauss widths from HERMES. At first glance, this seems to imply that the Gauss widths and EMC 
(where s = 525 GeV 2 ) are as large as at HERMES, while Fig. [TT] would suggest that they should be substantially 
larger. Indeed, if we extrapolate from Fig. [Til to EMC energies then we obtain for EMC about 2 times larger widths 
resulting in an about (30-40) % larger Cahn contribution to the cos ^-modulation. At this point, however, we have 
to recall that the Cahn effect description of this observable requires the neglect of quark-gluon-correlations in (|15l) . 
Such quark-gluon-correlations could be as large as (30-40) % compared to the contributions from quark correlators, 
as was found in [95[ in the context of a different observable. Thus, the intrinsic transverse momenta at EMC could 
be well compatible with the picture in Fig.[TTJ Future data from HERMES, COMPASS, and Jefferson Lab on A™^ 
will clarify the situation. 

To summarize: we conclude that the non-perturbative mechanisms responsible for intrinsic transverse parton mo- 
menta in DY and SIDIS are apparently compatible. On the basis of TMD-factorization and universality this is 
expected, but in view of the typically different energy scales probed in these processes it is not straight-forward to see. 



V. CONCLUSIONS 



In this work we discussed and reviewed the present understanding of intrinsic transverse parton momenta in SIDIS 
and DY. An important aim was to demonstrate that the popular Gauss model works very well in these processes. A 
similar in spirit study was presented in [6l| . But meanwhile many more data especially from SIDIS emerged, which 
allow more conclusive tests of the Gauss model, and make an update of previous results [6ll - l63j possible. 

More precisely, in the case of DY we have seen that the Gauss model is applicable if the transverse (dilepton) 
momenta much smaller than the hard scale. In SIDIS at energies probed at Jefferson Lab or HERMES we have seen 
that the Gauss model well describes all available data on cross sections or transv erse hadron momenta. It is interesting 
to remark that the Gauss model has received certain support from models [9(| Il29j and lattice QCD [l30j . 

We discussed the known fact that the Gauss widths increase with energy in DY. By comparing transverse hadron 
momenta at Jefferson Lab, HERMES and COMPASS we found indications for the energy dependence of the Gauss 
widths in SIDIS too. More precise data from SIDIS are needed, but on the basis of what is available now, we conclude 
that the intrinsic transverse parton momenta in SIDIS and DY are compatible. This demonstrates the universality of 
the pT-dependence of /f (x,pt) m SIDIS and DY and supports the factorization approach in terms of unintegrated 
correlators, although the support has presently a qualitative character. 

Our results are of importance for many practical applications. First, the Gauss model can now be used in SIDIS 
as an effective tool for the description of transverse parton momenta with more confidence. On the basis of presently 
available data we have seen no evidence for a worthwhile mentioning flavor- or x- or z-dependence of the Gauss 
widths. Second, we learned that the Gauss widths are energy dependent also in SIDIS, which is of importance when 
quantitatively comparing results from Jeff erson Lab, HERMES and COMPASS. The energy dependence of the Gauss 
widths in DY was well known before [1251 ]. Third, a solid understanding of pT-effects and their energy dependence is 
indispensable in order to study the azimuthal and spin asymmetries in DY and SIDIS. Our results will in particular 
be helpful to make estimates for the planned or proposed DY physics programs at COMPASS, GSI, J-PARC, U-70 
where the process will be probed at different energies. 

Especially in the context of the azimuthal asymmetries in unpolarized SIDIS a good quantitative understanding 
of intrinsic transverse momenta is of importance in the context of the Boer-Mulders effect. We have estimated that 
the associated cos(2<^))-modulation of the SIDIS cross section receives at Jefferson Lab, HERMES, and COMPASS 
energies sizeable 1/Q 2 power corrections from the Cahn effect — which is just one of the possible non-factorizing 
power corrections — in agreement with results from other studies [7l[ . Due to the larger Q in the available DY 
data the Cahn effect does not hamper the extraction of the Boer-Mulders function. 

As a byproduct, we confirmed the observation 62 1 that the Cahn effect can account for the EMC data [3l],[32j on the 
"twist-3" cos (^-modulation in the unpolarized SIDIS cross section, using the Gauss model parameters from HERMES. 
However, we were lead to the suspicion that this good description is likely to be due to a cancellation of effects due to 
the Gauss width broadening at the larger EMC energies and terms neglected in the so-called Wandzura-Wilczek-type 
approximation needed to justify here the "Cahn-effect-only" approximation. 
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To conclude, the Gauss model with carefully taken into account energy dependence of the Gauss widths is not 
only a convenient but also within a good accuracy well justified tool to describe intrinsic transverse parton momenta. 
This approach, especially now at the early state of art of studies of azimuthal and spin asymmetries in SIDIS and 
DY, represents a sufficient approximation for many practical purposes. The Gauss Ansatz remains to be tested, in 
particular, when polarization phenomena are included, and future data may demand to refine it, which will improve 
our understanding of intrinsic transverse parton momenta. 

Future steps, necessary when one will be interested in high precision and/or in going to high energies, for example 
DY at RHIC with <Js = 20 Ge V, will include the treatment within th e Collins -Soper-Sterman formalism [l2l ]. as 
implemented for instance in [l26| and consideration of scale dependence jl27l Il28| . 
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